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Abstract 

Entanglement of dipole-dipole interacting spins 1/2 is usually investigated when the energy of 
interaction with an external magnetic field (the Zeeman energy) is greater than the energy of dipole 
interactions by three orders. Under this condition only a non-equilibrium state of the spin system, 
realized by pulse radiofrequence irradiations, results in entanglement. 

The present paper deals with the opposite case: the dipolar interaction energy is the order of 
magnitude or even larger than the Zeeman one. It was shown that entanglement appears under 
the thermodynamic equilibrium conditions and the concurrence reaches the maximum when the 
external field is directed perpendicular to the vector connecting the nuclei. For this direction of 
the field and a system of two spins with the Hamiltonian accounting the realistic dipole-dipole 
interactions in low external magnetic field , the exact analytical expression for concurrence was 
also obtained. The condition of the entanglement appearance and the dependence of concurrence 
on the external magnetic field, temperature, and dipolar coupling constant were studied. 
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I. INTRODUCTION 



Appreciation of the role of quantum entanglement |1 3 1 as a resource in quantum telepor- 

n n n 

;ation [4|, quantum communication [5J, quantum computation [6(, and quantum metrology 
7|, |8|] has stimulated intensive qualitative and quantitative research. Entanglement, as the 
quantum correlation, can bring up richer possibilities in the various fields of modern technol- 
ogy. Therefore, in the past few years great efforts have been done to understand and create 
entanglement. Entanglement between two quantum systems can be generated due to their 
interaction only [l-3, 9]. It has recently been shown that, in a chain of nuclear spins s = 1/2, 
which is described by the idealized XY model for a spin system under the thermodynamic 
equilibrium conditions, entanglement appears at very low temperatures T w 0.5 /iK flo| . 

In most real quantum systems, such as dipolar coupling spin system, specific conditions 
for creation of the entangled states are requested. In two-and three-spin [ill and many- 
spin Q Clu8ter s of prol 8ubj e Ct ed t o Lon g magnetic ne.d, t— dipo.e-dipoL 
interactions and multiple pulse radiofrequence irradiations, the entangled state of a spin 
pair emerges at temperatures T « 20 mK. In these papers the cases were considered where 



the energy of interaction of the spins with the external magnetic fie 
is greater than the energy of dipole interactions by three orders 



d (the Zeeman energy) 
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121 ] . It was shown 



that at this condition only a non-equilibrium state of the spin system, realized by pulse 
radiofrequence irradiations, results in entanglement [l^, l^ . 



The present paper deals with the case opposite to those considered previously [Hi, |12j : the 
dipolar interaction energy is the order of magnitude or even greater than the Zeeman one. 
We investigate entanglement of two spins coupled by the realistic dipole-dipole interactions 
in a low external magnetic field under the thermodynamic equilibrium conditions. We study 
dependence of the critical temperature and magnetic field at which entanglement appears 
in this system on a dipolar coupling constant. 



II. HAMILTON! AN OF DIPOLAR COUPLING SPIN SYSTEM AND CONCUR- 
RENCE BETWEEN NUCLEAR SPINS 1/2 

Let us consider a system of N spins coupled by long-range dipolar interactions and 
subjected to an external magnetic field, H = H z. The total Hamiltonian of this interacting 
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system can be written as 

H = H Z + H dd (1) 

where the Hamiltonian H z describes the Zeeman interaction between the nuclear spins and 
external magnetic field (here we used H — 1) 

N 



H z =u J2lk, (2) 



ujq = ■jHq is the energy difference between the excited and ground states of an isolated spin, 
7 is the gyromagnetic ratio of a spin, I| is the projection of the angular spin momentum 
operator on the z- axes. The Hamiltonian Hdd describing dipolar interactions in an external 
magnetic field jl^ : 
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H dd = ^2 ~{ i 1 ~ 3 cos2 ^' fc ) 



| sin 26 jk [e-*** (///+ + Ip z k ) + e*» (Iff + - \ sin 2 9 jk [e'^Ipt + e 2 ^ljl k ] } 

(3) 

where rj k , 6j k , and ipj k are the spherical coordinates of the vector fj k connecting the j—th 
and k— th nuclei in a coordinate system with the z-axis along the external magnetic field, 
H , If and J^" are the raising and lowering spin angular momentum operators of the j-th 
spin. We consider the situation when it is necessary to take into account all the terms of 
the Hamiltonian of the dipole-dipole interactions, and not trusnckete any ones. 

In the thermodynamic equilibrium the considered system is described by the density 
matrix 

P = Z-'e W (-JL), (4) 



where Z = Tr {exp (—H/hsT)} is the partition function, is the Boltzamnn constant, 
and T is the temperature. We will analyze entanglement in the spin system described by 
the density matrix (4). 

n 

In order to quantify entanglement, the concurrence C is usually used [15]. For the maxi- 
mally entangled states, the concurrence is C — 1, while for the separable states C — 0. The 
concurrence between the quantum states of two spins presented in the Hilbert space as a 



matrix 4 x 4 is expressed by the formula 
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C = max I 0, 2A - ^ A< \ (5) 



i=l 




where A = max {Ai, A2, A3, A4} and Aj (i = 1, 2, 3, 4) are the square roots of the eigenvalues 
of the product 

R = pp (6) 

with 

P = (<Ty <8> <Ty) P {<Jy ® <T W ) (7) 

where p the complex conjugation of the density matrix (4) and a y is the Pauli matrix 



(8) 



III. ENTANGLEMENT IN PAIR OF SPINS 

We examine dependence of the concurrence, C, between states of the two spins 1/2 on the 
magnetic field strength and its direction, dipolar coupling constant, and temperature. The 
numerical calculation of entanglement of the spins at arbitrary orientation of the magnetic 
field are performed using the software based on the Mathematica package. The results of the 
numerical calculation show that concurrence reaches its maximum at the case of 9 — f and 
(p = (Fig. 1) and we will consider this case below. This orientation of the spins allows us 
to obtain for concurrence as an exact analytical function of the temperature, magnetic field 
and dipolar coupling constant 7 2 '/rf 2 - Using the exact diagonalization of the density matrix 
(4) we obtain the concurrence in the following form: 

A: — A_ — ei cosh - A 
C <!) max { 0, — -± — * } , (9) 

cosh I + cosh ^ — - A 



where 



16g + 9<P cosh ± Qd^^^^lQP + 9^ cosh 2 

\ 16/3 2 + 9d 2 



(10) 



with 13 = and d = ^r^. 

At high temperature and low magnetic field (j3 « 1) and/or small dipolar coupling 
constant (d « 1) the expression in the figure brackets (9) becomes negative and, therefore, 
entanglement is zero. Equating this expression to zero we obtain the critical parameters: 



temperature T c , strength of magnetic field H C} and dipolar coupling constant at which the 
entanglement appears in a spin pair. Figure 2 presents the phase diagram in which the 
boundary between the entangled and separated states is determined from equation: 



A + - A_ - ei cosh ^ = 0. (11) 

For example, at d = 1 entanglement can be achieved at /3 > 2.26. To estimate the critical 
temperature let us consider fluorine with 7 = 4. 0025^P and the dipolar interaction energy 
typically of order of a few kHz (in frequency units) |l4j. Taking H = 3 G we have cuq = 12 
kHz, which leads to T c = 0.33 fiK. The estimated value of temperature is in good agreement 
with those reported early for the spin system s — 1/2 with the XY Hamiltonian in absence of 
a magnetic field [IQJ]. Both models give also qualitatively similar dependences of concurrence 
on temperature but the model considered by us predicts appearance of entanglement in an 
external magnetic field higher than the critical value. It is interesting that the ordered states, 
such as antiferromagnetic, of nuclear spins were observed in a calcium-fluoride CaFo single 
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crystal at T = 0.34 /iK |14J . Il6j . This structure is characterized by domains in the form of 
layers perpendicular to the external magnetic field. The magnetization inside the crystal 
is parallel to the external field and reverses its sign from a layer to a layer, while the total 
magnetization is zero. It is well know that when the magnetization reaches the maximum, 
all the spins are aligned along the field and entanglement is absent 

as. 

Entanglement 

can appear if magnetization is less than its maximum. Therefore, it is reasonable to assume 
that simultaneously with the transition to the ordered state there arises entanglement of 
spins from different layers. 

Figure 3 shows concurrence as a function of both parameters /3 and d at 9 — | and 
ip = 0. At large temperature and low magnetic field concurrence is zero. One can see that 
the concurrence increases with the magnetic field and inverse temperature and reaches its 
maximum. Then the concurrence decreases. Figure 4 shows the concurrence as a function 
of the magnetic field at a constant temperature, (Fig. 4a) and as a function of the inverse 
temperature at a constant magnetic field (Fig. 4b). In the both cases concurrence remains 
zero up to a certain value of the magnetic field (Fig. 4a) or of the inverse temperature (Fig. 
4b), which depends on the coupling constant. The following increase of the magnetic field 
or inverse temperature leads to eventually rising. In the case of increasing the magnetic field 
concurrence increases up to the maximum and then decreases as a magnetic field increases 



(Fig.4a). Another behavior is observed at an increase of the inverse temperature (Fig. 4b): 
concurrence monotone grows with 1/T up to a steady value depending on a magnetic 
field and the following increase of the inverse temperature does not cause any change of 
concurrence. 



IV. DISCUSSION AND CONCLUSIONS 

It was obtained that in zero magnetic field the system is in a separable state. The 
system becomes entangled when the interaction energy of spins with the magnetic field are 
of the order of the dipolar interaction energy. Then, with increasing magnetic field the spin 
state tends to separable one. At a small dipolar coupling constant (d<l) from the exact 
analytical solution (9) we obtain the following expression for the concurrency 



C = max < 0, 



1 



2 cosh 2 § 



(12) 



Therefore, at these conditions the states of the system are always separable, C — 0. Entan- 
glement appears in the course of increasing the dipolar coupling constant. To distinguish an 
entangled state from separable ones, it is important to determine an entanglement witness 



(EW) applicable to the considered quantum system 
one of the main problems of the experimental study o 
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18| . The determination of EW is 
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19], magnetic susceptibility [20[, magnetization [21 
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the entangled states. Internal energy 



22] , and intensity of MQ coherences 



24j were used as EW in different quantum systems. With the aim to obtain the 



correlation between the nuclear magnetization and concurrence, using (4) and the definition 
of nuclear magnetization M z = Tr (pl z ), we obtain the exact expression for magnetization 
as a function of parameters and d at 9 12 = | and tp 12 = 0: 



M r 



-4/3sinhV^±^ 



(16/3 2 + 9d 2 ) (^cosh v^/W + ef cosh \ 



(13) 



As example, at d — 3, the relation between the concurrence and magnetization can be fitted 
by C = -0.71(M + 0.26) (Figure 5). 

Concurrence, the measure of entanglement between the states of the two spins, depends 
on the orientation of the magnetic field relative to vector r connecting the nuclei. At 9 = 



and 7r the states are separable and the concurrence reaches its maximum at the case of 9 — | 
(Fig. 1). This effect can open a way to manipulate with the spin state by a rotation of the 
magnetic field or a sample. 

In conclusion, investigation of entanglement in a spin 1/2 system under the thermody- 
namic equilibrium conditions showed that the entangled state can be achieved by application 
of a low external field when the Zeeman interaction energy is the order of or even less than 
the dipolar interaction one. It was estimated that for magnetic field H = 3 G, the entan- 
gled state in a two-spin system arises at temperature T < 0.33 /iK. The correlation between 
concurrence and nuclear magnetization is considered and it was shown that concurrence is 
well fitted by a linear dependence on the magnetization in the temperature and magnetic 
field range up to a deviation of the magnetization from Curie's law (j3 = 3.32, Fig. 5) . 
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A. Figure Captions: 

Fig. 1 (Color online) Concurrence as a function of the parameter /3 = uo/ksT and 
magnetic field direction at ip = and d = 3. 

Fig. 2 The phase diagram. Line presents boundary between the entangled and separated 
states in the plane (3 — d. 

Fig. 3 (Color online) Concurrence as a function of the ratios of the magnetic field strength 

2 

(a>o) and dipolar coupling constant (^-) to ksT. 

r 12 

Fig. 4 (Color online) Concurrence vs. magnetic field at T = const (a) and vs. tempera- 
ture at H — const (b) for various dipole interaction constants, (a): black solid line - d — 0.5; 
red dashed line - d = 2 ; blue doted line - d = 10. Magnetic field is given in units of . 
(b) black solid line -4 = 3; red dashed line -4 = 5; blue doted line - 4 = 10. Temperature 

2 

is given in units of 3 7 , . 

Fig. 5 (Color online) Absolute value of magnetization (black solid line) and concurrence 
(red dash line) as a function of (3 = Fitting of the concurrence (blue dash-dot line) by 
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C = -0.71(M + 0.26) at d = 3. 
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